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ABSTRACT 


The Institute of Industnal Science, Tokyo, Japan, initiated the Chiba downhole array, in 1982. 
This was instrumented with several densely spaced free field and downliole accelerometers, thus 
constituting a three dimensional array. In this study, the seismic site response during 27 
earthquakes recorded at the site with varying peak ground accelerations, frequency content and 
duration is analyzed. Correlation analyses of the recorded downhole accelerations are performed 
to calculate the average shear and primary wave velocity profile. A two-layer solution was also 
employed to calculate the average shear wave profile. In addition, the seismic shear stress-strain 
histories are directly calculated from the downhole lateral accelerations. These histories are then 
used to estimate the soil stiffiiess and material damping in the linear, small strain range. The 
identified shear wave as well as primary wave velocity profile was compared witli the 
geophysical measurements and with those reported by other researchers. The observed non- 
softening of the soil below 10 m depth even during the strongest recorded earthquake was 
explained. 
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Chapter 1 

INTRODUCTION AND OBJECTIVES 

1.1 Introduction 

Structural damage due to seismic excitation is directly related to the underlying ground 
response characteristics. Site condition, geology, soil types and soil pore water interaction 
affect these characteristics. Recent major earthquakes, such as those at Mexico (Sept. 19, 
1985), Loma Prieta (Oct.l7, 1989) and BhuJ (Jan.26, 2001) emphasize the effects of site 
condition on the resulting ground surface response, and consequently on the extent of 
seismically induced damage. Several earthquake resistant design and analyses problems 
require dynamic / cyclic soil properties and site characterization including, 

• Local site response (site amplification) 

• Seismic microzonation for hazard mitigation 

• Liquefaction and densification of level or almost level sites 

• Soil structure interaction 

• Geotechnical structure response 

• Slope stability 

• Site improvement. 

All these problems require knowledge of cyclic soil properties defined and measured along 
characteristic loading paths of interest. Site amplification / local site response (wave 
propagation problem) is important by itself and is also the starting point for the other 
problems listed above. Soil amplification characteristics can provide information for the 
prediction of ground motion for a specific site and for the seismic safety assessment of 
existing structures. 

Hence, investigation of soil amplification characteristics and estimation of dynamic soil 
properties of sub-surface soil is extremely important as local site conditions can significantly 
influence the characteristics (amplitude, frequency content, and duration) of earthquake 
ground motion and the degree of damage caused by it. 
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Assessment of site amplification effects necessitates a thorough understanding of the 
associated soil dynamic response characteristics and mechanisms. In-situ testing procedures 
currently followed, provide convenient means of measuring soil low-strain properties such as 
shear wave propagation velocities. At larger strain levels, soil properties are often evaluated 
through laboratory testing. However, due to limitations in undisturbed soil sampling 
techniques and difficulties in reproducing in-situ stress-strain and seismic loading history, 
data provided by laboratory tests may not reflect the site actual behavior, especially when 
subjected to strong-motion excitations. 

Recorded site seisimc response is a valuable complement to in-situ and laboratory 
geotechnical investigation techniques. In particular, analysis of downhole seismic array data 
provides unique information of actual soil (and overall site) behavior over a range of loading 
conditions that are not readily covered by in-situ or laboratory experimentation procedures. In 
view of the need for a better understanding of site dynamic response, downhole 
instrumentation is being increasingly deployed at seismically active locations especially in 
the United States, Japan, Mexico, and Taiwan. The recent fast accumulation of earthquake 
records firom seismic arrays in well-documented sites provide a powerful new tool to measure 
and investigate the soil properties of interest (site characterization / amplification 
characteristics) under actual earthquake condition. 

An earthquake may be viewed as a full-scale experiment and recorded input and output may 
be used to identify the stiffhess and damping properties, and is a valuable complement to 
other geotechnical investigations. Evaluation of soil properties from earthquake records 
provide information on soil behavior over a range of excitation that are not readily covered 
by in-situ or laboratory experimentation and more importantly without resorting to in-situ 
sampling (Kramer et al. [32]), which may not be a very representation of the overall 
composition. 

Once the array station established and seismometers installed, sequential earthquake records 
can be obtained for further analysis, unlike as in the case of laboratory analysis where the 
samples would be required every time from the required depth at which the in-situ test results 
are needed. Further the analysis done on borehole recordings may be correlated with results 
of laboratory testing. 
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In 1 982, a unique three dimensional array system was installed at Chiba experiment station of 
Institution of Industrial Science, University of Tokyo. A database of acceleration time history 
of 27 earthquake events was made available for various boreholes (CO, Cl, C2, C3, C4, PI, 
P2, P3, P4, P5, P6, P7, P8, P9, and PO). The site location was selected because of its intense 
seismic activity and the presence of thick soil strata above bedrock. In the present study, 
acceleration time history records of the borehole CO, where the Peak Ground Acceleration 
(PGA) is greater than 10 cm/sec^, are used for estimating the dynamic soil properties. These 
earthquakes were found to vary in peak acceleration, shaking duration and frequency content 
bandwidth; thus providing a valuable source over a broad range of excitations. 

1.2 Objectives 

The main objectives of this study were to analyze the Chiba seismic records and to identify 
the site dynamic response properties and mechanisms. Free field surface and downhole 
acceleration records were used, 

• to identify the site resonant frequencies 

• to evaluate the site shear wave propagation characteristics 

• to document the presence of nonlinear response at large Peak Ground 
Acceleration (PGA) levels and 

• to evaluate the variations of soil stiffiiess and damping properties. 

1.3 Methodology 

Identification of site dynamic response properties and mechanisms from seismic records 
requires preliminary analyses to assess data quality and classify the recorded response and 
loading conditions. The identification methodology consists of: 

• Pattern recognition and classification of the recorded response. 

• Selection of a class of appropriate analytical / computational models and 
assessment of the most influential djmamic response parameters. 

• Estimation of these parameters. 
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In 1982, a unique three dimensional array system was installed at Chiba experiment station of 
Institution of Industrial Science, University of Tokyo. A database of acceleration time history 
of 27 earthquake events was made available for vanous boreholes (CO, Cl, C2, C3, C4, PI, 
P2, P3, P4, P5, P6, P7, P8, P9, and PO). The site location was selected because of its intense 
seismic activity and the presence of thick soil strata above bedrock. In the present study, 
acceleration time history records of the borehole CO, where the Peak Ground Acceleration 
(PGA) is greater than 10 cm/sec^, are used for estimating the dynamic soil properties. These 
earthquakes were found to vary in peak acceleration, shaking duration and frequency content 
bandwidth; thus providing a valuable source over a broad range of excitations. 

1.2 Objectives 

The main objectives of this study were to analyze the Chiba seismic records and to identify 
the site dynamic response properties and mechanisms. Free field surface and downhole 
acceleration records were used, 

• to identify the site resonant frequencies 

• to evaluate the site shear wave propagation characteristics 

• to document the presence of nonlinear response at large Peak Ground 
Acceleration (PGA) levels and 

• to evaluate the variations of soil stiffness and damping properties. 

1.3 Methodology 

Identification of site d 3 mamic response properties and mechanisms from seismic records 
requires preliminary analyses to assess data quality and classify the recorded response and 
loading conditions. The identification methodology consists of: 

• Pattern recognition and classification of the recorded response. 
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• Estimation of these parameters. 
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The conducted analyses fall within the pattern recognition phase. They are based on non- 
parametiic identifications to classify the recorded seismic response at Chiba and evaluate the 
most influential response parameters and mechanisms. 


1.4 Analyses Techniques 

As stated above, the objectives of this study were to analyze the recorded seismic soil 
response at Chiba array and to identify the site dynamic response properties. Two classes of 
techniques were employed to attain these objectives: 

• Correlation and spectral techniques to estimate site low-strain parameters such as 
shear wave velocities, resonant firequencies and to assess the effects of non- 
linearities on these parameters. 

• A simple non-parametric identification technique was developed to estimate the 
site shear stress-strain histories and to evaluate variations of stiffiiess and 
damping. 


Spectral analysis was based on the Fourier transforms of the recorded data (using a fast 
Fourier transform algorithm). Welch technique was employed in estimating energy spectra to 
reduce random errors. A new non-parametric identification procedure recently developed was 
used to evaluate seismic shear stress-strain histories directly firom downhole acceleration 
records. 

This procedure is simple to implement and is well suited to site geometry and downhole 
accelerometer spacing. Furthermore, it does not require the availability of bedrock (input) 
motions, in contrast to many conventional techniques. At Chiba site no input records are 
available. In estimating stress and strain histories, filtering was performed using non- 
recursive zero-phase digital time-domain filters through Harming window. 
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1.5 Organization of the Report 


In this study report, site geotechnical profile geophysical measurements at borehole CO and 
acceleration data of Chiba array are presented in Chapter 2. Literature review covenng 
similar research work in the area of downhole arrays is presented in chapter 3. Shear and 
primary wave velocities estimated using a non-parametnc identification technique (cross- 
correlation method) are compared with geophysical measurements in chapter 4. 

Recorded acceleration time histories are used to identify the stiffness and damping properties 
of in-situ soil under actual earthquake conditions and the computed shear stress-strain 
histories are presented in chapter 5. Using a two-layer solution (a non-parametric 
identification), shear wave velocity profile between two downhole stations is estimated and 
compared with geophysical measurement profile in chapter 6. Summary and conclusions 
based on the results of three analyses tecliniques are presented in chapter 7. 
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Chapter 2 

STRONG MOTION DATA - CHIBA ARRAY 


2.1 Introduction 

One of the state-of-the-art technologies in strong ground motion measurement is array 
installation and monitoring the data. The array observation system is an ideal tool for 
accumulating the strong-motion data. An array system provides information regarding strong- 
motion parameters viz. acceleration, velocity and sometimes strain by using systematically 
placed instruments oriented deep below the ground surface to monitor the peak earthquake 
signals as per the requirement. 

The aim of observation dictates the arrangement of seismometers. One of the typical types is 
a two-dimensional surface array SMART- 1 on thick alluvial site at Lotung, Taiwan and 
SMART-2 array on hard site in the Haulien area (Yeon-Tein et al. [51]). Another 
arrangement used for the present analysis is a dense seismograph array network installed at 
Chiba experiment station (Nagata et al. [36]). In this eirray, seismometers are placed very 
densely both on the ground surface and in the boreholes, constituting a three-dimensional 
network. The records obtained from this array are useful in evaluating the spatial correlations 
and in examining the amplification theory. 

Arrays are also useful for studying earthquake response mechanisms and spatial distribution 
of many important earthquake parameters such as, 

• Amplitude parameters (Peak acceleration, peak velocity, peak displacement. 
Sustained maximum acceleration and velocity, Effective design acceleration) 

• Frequency content parameters (Ground motion spectra, Fourier spectra. Power 
spectra, and Response spectra) 

• Spectral parameters (Predominant period, Bandwidth, Central fi-equency. Shape 
factor) 

• Duration. 



In addition, these records offer a solid basis for, 


• Calculation of in-situ dynamic properties 

• Refinement of empirical and computational predictive techniques 

• Calibration of numerical models 

The Chiba array has been successfully operating since 1982 and over 160 events have been 
recorded. A database has been created for the effective utilization of these earthquake 
records, which can be obtained from University of Tokyo, Japan. 

2.2 Geotechnical profile 

The Chiba experiment station is located about 30 km east of Tokyo as shown in Figure 2.1. 
The longitude of the station is 140° 6' 37" E and the latitude is 35° 37’ 17" N . 



W 
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Figure 2.1: Location of Chiba experiment station 




The layout of the boreholes instrumented with accelerometers is shown in Figure 2.2. The co- 
ordinates of each borehole at ground surface are listed in Table 2.1. The data obtained from 
borehole CO is used in the analysis for the present work. Seismometers are placed in the 
boreholes at various depths and the location of borehole accelerometers is given in Table 2.2. 
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Figure 2.2: Layout of boreholes 


Table 2.1: Co-ordinates of boreholes at ground surface 


Borehole 

Relative to CO 

Elevation (m) 

East (m) 

East (m) 


0.00 


13.52 

Cl 

3.54 

3.54 

13.52 

C2 

-3.54 

3.54 

13.52 

C3 

-3.54 

-3.54 

13.52 

C4 

3.54 

-3.54 

13.52 

PI 

15.00 


13.52 

P2 

0.00 

15.00 

13.52 

P3 

-15.00 

0.00 

13.52 

P4 

0.00 

-15.00 

13.52 

P5 

-45.13 

-115.38 

12.59 

P6 

-128.92 

-61.42 

12.64 

P7 


4.67 

13.65 

P8 

-216 13 


11.91 

P9 

-109.01 

178.74 

. 10.92 

PO 

44.94 


12.10 
























Table 2.2: Location of borehole accelerometers 


Depth 
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Extensive geotechnical site characterization was carried out at these locations in terms SPT 
and geophysical measurements. Topographical and geological conditions of the site are 
generally simple, as the ground surface is almost flat. Figure 2.3 shows the soil profile of 
borehole CO and the elastic wave velocities measured by downhole shooting method. 



Figure2.3: Soil profile and geophysical measurements at borehole CO 























Top 3 - 5 m of the site contained loam with standard penetration (N) value less than 10 and 
sandy clay underlies this loam layer with a thickness of 2 - 4 m, again with N values less than 
10. A diluvium sand layer lies under the clay layer with N values greater than 20-30. The 
sand layer (whose stiffness increased with depth) is interspersed with clayey layers of 
relatively small N values. In spite of slight differences in the depth of boundaries between 
different layers from one borehole to another, an overall agreement indicates a relatively 
simple soil structure. The ground water table was found to be lower than GL-5.0 m. 


2.3 Instrumentation 

A piezoelectric t 5 ^e acceleration transducer was used in the array observation. Three 
transducers (two horizontal in East-West (EW), North-South (NS) and one in Up-Down (UD) 
directions) and their amplifiers are installed in cylindrical steel casing with an external 
diameter of 65 mm and a length of 335 mm. The accelerometers were installed in boreholes 
with diameters of 116 mm. Maximum of five accelerometers in each borehole was installed 
and fixed at predetermined depths (1, 5, 10, 20 & 40 m) using cement mortar. 

Signals from the seismometers (placed in the boreholes to record signals in three directions 
viz. EW, NS and UD) were recorded at every 0.005 seconds by three units of 64 channel 
digital recorder. The system is always kept in a full operational status. Recording devices are 
activated the moment any one of the three components at p540 (GL - 40 m in p5 borehole) 
exceeds an acceleration of l.Ocmsec"' . 

The system continues in operation for 30 seconds after the motion falls below this trigger 
threshold level. The recorder consists of a magnetic tape with a recording capacity of 30 
minutes. Timing information is internally generated, and in addition the absolute time is 
corrected hourly by using the signal from NHK radio station. Detailed layout of 
accelerometers at different depths in the borehole CO is shown Figure.2.4. 
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® Im depth accelerometer 
® 5 m depth accelerometer 

® 10 m depth accelerometer 

@ 20 m depth accelerometer 

@ 40 m depth accelerometer 


Figure 2.4; Layout of accelerometers at borehole CO 


2.4 Database development of Chiba array 

The Chiba array database consists of the data on major earthquakes whose peak ground 
acceleration at CO (GL — 1 m) is larger than 20c7nsec"^ or whose maximum steel pipe strains 
are larger than 5x10“®. The basic information on included earthquake records (only 25 events 
are considered) of the Chiba array database is listed in Table 2.3. Voltage signals from the 
seismometers are stored in digital magnetic tapes. These voltage values are multiplied with 
the conversion factor to give the values of acceleration, velocity and strain. Baseline 
correction is applied in this stage to have zero-mean for acceleration records. 



If a seismograph does not start until a triggering level of motion is reached, the entire 
accelerogram is in error by the level of motion at the time of triggering. Integration of an 
uncorrected acceleration time history will produce a linear error in velocity and quadratic 
error in displacement. Comparing the same signals recorded in the tliree tape units, the time 
lag correction is applied. The recording time is sometimes too long to store a whole record in 
the database. Then the recording duration is truncated at the time when the cumulative power 
of horizontal acceleration records at CO exceeds 99% of the total recording time. After these 
operations, the strain and velocity records are stored in the database. However, further 
corrections are applied to the acceleration records and a high cut-off filter is employed 
considering the sensitivity of the accelerometer. 


Table 2.3: Basic information on earthquake records 
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The installation of seismographs in the borehole is very important because they are installed 
at great depths (up to 40 meters) in the boreholes without any visual observation. While 
setting the seismograph in borehole, utmost care is taken to give it the correct orientation, but 
still orientation errors do occur and were reflected in the seismograms. To correct these 
erroneous seismograms, orientation errors for the faultily oriented instruments are estimated 
and the records are corrected. 

2.5 Organization of Database 


Two databases, one for acceleration and the other for strain and velocity are stored in 
magnetic tapes. Each database consists of 27 sequential data sets corresponding to the 27 
events as shown in Figure 2.5. Each data set includes 135 records for the acceleration 
database and 37 records for the strain and velocity database. 
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Figure: 2.5 Organization of Chiba array database 









Each time history data is composed with a heading data that includes the event code (EEQK), 
the component code (ICMP), the title of the record, the maximum value of the record, the 
time interval (sampling interval) and the number of time steps. The details of the event and 
component codes are shown in Figure 2.6 and the combination code is unique for each 
record. The time history data is followed after the heading data. 


lEQK 

1 coL 

o 

o 

CM 1 

3 coL 

4 col 


year event number 
in each year 


ICMP 


• 





o 

o 

O 

O 

o 

u 

o 

u 

u 

o 


CM 



ITS 

; 

: 

: 


t 

! 

L. 


1 :acceleration 
2: velocity 
(3: displacement) 


4: ground strain 
5:pipe strain 

6: relative displacement of pipe joint 
7; time code 


1 : NS comp . 
- 2 : EW comp . 
3 : UD comp . 


4:building axi$~y 
5:building axls-x 
6:right side of pipe 


7:left side of pipe 
8:upper face of pipe 
9: lower face of pipe 


0: GL - Im 
1: GL - 5m 
2: GL “10m 
4; GL “20m 
8: GL “40m 


10: borehole CO, 
llrborehole Cl. 
12: borehole C2, 
13:borehole C3, 
14:borehole C4, 
21: borehole PI, 
22: borehole P2, 


23: borehole P3 , 
24:borehole P4 , 
25: borehole P5 . 
26:borehole P6 . 
27: borehole P7 , 
28: borehole P8 . 
29 ‘.borehole P9 , 


30:borehole PO 
40: floor acc. 

50: floor vel. (*10) 

55 : floor vel - ( *1 ) 

60“B9 : ductile iron pipe 
70“79:steel pipe 
80-89 ‘.ground 


Figure 2.6: Event and Component codes 
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2.6 Data Procurement / Retrieval 


Data records of 27 events recorded at various boreholes are made available in a magnetic 
tape storage device in the compressed form (Zip format). The records at 1,5, 10, 20 & 40 m 
depths in EW, NS and UD directions at borehole CO were separated from each data set. Using 
SUN Solaris SPARC work station at IIT Delhi, each compressed data set of the database is 
unzipped and the required data records are copied. 5-10 Mb hard disc space is required for 
each data set (after unzipping) and is about 80 Mb for all the records, hence sufficient 
memory space should be ensured. Storage and management of data records is very essential 
before proceeding to the analysis work when the disc space is constrained and CD-ROMs can 
come in handy. 

2.7 Plots of Selected Database and Observations 

Five different events were selected for plotting the acceleration time histories and Fourier 
spectra as shown in Figures 2.7 - 2.16. These events are: 

• 8706 where PGA is 14 cm/sec^ in NS direction 

• 8903 where PGA is 28.9 cm/sec^ in NS direction 

• 8519 where PGA is 82.2 cm/sec^ in NS direction 

• 8525 where PGA is 75.7 cm/sec^ in EW direction 

• 8722 where PGA is 327.1 cm/sec^ in NS direction 

The amplitude of acceleration is observed to be higher at surface than at 40 m depth as can be 
seen from the Figure 2.7-2.16. The acceleration at the ground surface is about 3-4 times 
more than the acceleration at the 40 m depth. The Fourier spectrum shows that the frequency 
content of the record at 40 m depth and at the surface are different and also the energy of 
different waves is being considerable only up to a maximum frequency of 10 Hz. The energy 
and frequency content of the wave is getting modified gradually as it reaches the surface. 
Therefore it can be said that the soil behavior under earthquake loading at the ground surface 
is different than the soil at other depths, also considerable changes in acceleration and 
frequency content could be observed mostly in surface and 5 m records. 
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Figure 2.7: Acceleration time histories during the earthquake event 8706 
at borehole CO 
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Figure 2.11: Acceleration time histories during the earthquake event 8519 
at borehole CO 
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Figure 2,13: Acceleration time histories during the earthquake event 8525 
at borehole CO 
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Figure 2.16: Fourier Spectra during the earthquake event 8722 at borehole CO 









Chapter 3 

LITERATURE REVIEW 


3.1 Introduction 

Parameter estimation of soils [1, 2, 18, 33, 34, 45, 52, 53, 54, 55] and structures [6, 13, 16, 
19, 31, 37, 39, 41] from dynamic response data has enjoyed considerable attention in recent 
years. The techniques that have been developed range from linear deterministic models to 
nonlinear stochastic models. The applications vary from improving models for damage 
detection to identifying inputs to a system in order to control its response. This notable surge 
in related research activity is largely due to recent trends towards increased availability of 
computational resources and availability of experimental data and earthquake records. 

In this chapter, a brief bibliographic review is outlined. This review is intended to give an 
overall idea about the research activities in the fields that are directly related to this 
investigation i.e., identification of dynamic soil properties using downhole array data. Work 
related to estimation of dynamic soil properties of an earth dam using strong motion data, 
fimt of this kind, is presented in section 3.2. 

Several efforts in recent years to use downhole array data for parameter estimation are 
compiled in section 3.3. Research work that was reported on Chiba array data is compiled in 
section 3.4. Section 3.5 describes other studies conducted on downhole seismic data, that 
may not directly be related the current work, for the sake of completeness. A brief summary 
is presented in section 3.6. 

3.2 The First Attempt 

Research efforts to identify soil behavior from earthquake records are scarce. An early 
success at estimating dynamic soil properties using measured earth dam response was 
reported by Abdel - Ghaffar and Scott [1,2]. ' 



A shear beam model was used to estimate the shear moduli and damping ratio for the gravely 
and rolled-fill Santa Felicia earth dam. The soil was modeled as a hysteretic SDOF system 
with a nonlinear restoring force, and narrow band-pass filtering of the data isolated 
fundamental frequency of dam. An estimate of shear modulus and damping (Figure. 3.1) was 
obtained by plotting absolute accelerating of the dam against crest displacement relative to 
the base. 



Figure 3.1: Hysteretic stress-strain relationships used to calculate the shear moduli 
and damping factors of earth dam (Abdel - Ghaffar et al. [2]) 


Recently, Lin and Chao et al. [33] suggested guidelines to improve the hysteresis loops, that 
were observed to be distorted in Abdel - Chaffer et al. [2] study, based on a numerical 
simulation on earth dam as a non-linear, purely hysteretic SDOF sfructure subjected to a 
given base motion. From the calculation of shear moduli and damping factors, it was 
observed that the loops obtained were undistorted for the same time window in comparison 
between Abdel - Chaffer analysis and present analysis. It was concluded that to minimize 
errors, only those loops of larger strain amplitudes only should be used. 




3.3 Work related to Downhole arrays 


Strong motion earthquake records provide one of the few sources of information on the 
dynamic behavior of soil / overall site behavior at large amplitude response. Under 
earthquake conditions, soil may sustain unpredicted damage or may display unanticipated 
strengths; such information is not available through other means. This has motivated an 
extensive deployment of instrumentation to record soil response / site amplification in 
seismic regions. Currently, high quality records of well-instrumented soil sites (Japan, USA, 
Mexico, and Taiwan) are available, which are being used for a verity of research 
observations. This section describes some of these efforts related to current study. 

Zeghal and Elgamal [53] used a non-parametric identification technique to evaluate the 
mechanism of wildlife refuge site liquefaction. Surface and downhole accelerations and pore 
water pressures recorded during the 1987 earthquakes were used to obtain direct estimates of 
average shear stress-strain and effective stress path histories, which provided valuable insight 
in to the site seismic behavior during liquefaction and associated loss of stiffiiess. 

Elgamal et al. [14] used the cross-correlation method to calculate the shear wave velocity and 
propagation mechanism using Lotung free-field surface and downhole acceleration records. 
Resonant and modal configurations were estimated using spectral analyses. It was concluded 
that spectral and correlation techniques may be utilized to identify soil low-strain properties 
and site resonant characteristics. 

Matsunaga et al. [35] adopted cross-correlation method and ray theory for calculating time 
lag and phase velocities of body waves using dovrahole array records established in the 
Tokyo international airport, Japan. In both the studies, the computed S wave propagation 
velocity values are in good agreement with the measured values (geophysical tests). 

Campanella et al. [8] has given the frequency domain calculation to obtain the timelag for the 
wave to travel between the two adjacent stations. Crossover methods and Cross-correlation 
methods were applied on signals from downhole seismic cone analysis. It- was established 
that the Cross-correlation method of obtaining the time lag is more reliable than the 
Crossover method. 



Huang and Chiu [21] analyzed four-level Dahan downhole array recordings (at 0, 50, 100 and 
200 m depth) to investigate site amplification due to the near-surface structure. They used 
eight sets of data to study site amplification. Five techniques of spectral ratios have been 
tested to determine the most reliable way to estimate site response. The study reported the 
superiority of Ratio of Root-Mean-Square (RMS) Fourier amplitude spectra. Results show 
that the amplification factors are nearly uniform at the depths of 50 and 100 m and agree well 
with those from time-domain analyses and the transmission-coefficient estimation. 

Bersenev et al. [7] studied the soil nonlinear behavior under dynamic loading conditions and 
carried out spectral ratio analysis on weak and strong motion data at SMARTl and SMART2 
array at Lotung and Haulien. In SMARTl analysis, amplification of soft alluvial site relative 
to the rock and in SMART2 analysis, amplification of terrace deposits with respect to 
alluvium deposits was calculated. In both the cases, significant differences between the weak 
and strong motion application functions were observed. When the thickness of the surface 
layers is significantly less, the low frequency contribution from nonlinearity was less 
important and long period waves were not observed. 

Honjo et al. [20] given Extended Bayesian Method (EBM) to study the propagation 
characteristics of earthquake motion from bedrock. The accuracy of the method was 
compared with the conventional Least Square Method (LSM). The array records available 
from downtown Tokyo, Japan were used in the analysis. Transfer functions were smoothened 
using Parzen window. The estimated shear modulus and damping ratio by inverse analysis 
(EBM) were in good agreement with the observed records, but not in the case of conventional 
analysis (LSM) and hence proved the effectiveness of the EBM over LSM. 

Chang et al. [11] estimated the in-situ dynamic properties to provide field evidence of non- 
linear soil behavior during earthquake excitation and to evaluate the accuracy of dynamic 
properties obtained firom geophysical measurements and laboratory tests. From the Lotung 
downhole array data, shear wave velocities were estimated firom the fundamental resonant 
frequencies identifiable from Fourier spectral ratios. Non-linear behavior of shear strain was 
observed, which was estimated by linear ground response deconvolution analysis based on 
inferred shear wave velocity profiles. Amalysis on low excitation and high excitation levels 
were calculated in different time windows. The back-calculated shear wave velocity, moduli 



values were compared with the velocity obtained from geophysical tests and shear modulus 
reduction curves respectively. Substantial drop in shear wave velocity at high excitation 
levels was observed, when compared with the shear wave velocity from geophysical tests and 
hence the results clearly indicated the non-linear behavior of soil during strong motion 
earthquake. 

Directly assessing the non-stationaiy data, Glaser [18] found that parametrie modeling of the 
system is superior to the traditional Fourier methods, to estimate in-situ soil properties, 
especially when soil is subjected to large strains. 

3.4 Studies on Chiba Downhole array 

Katayama et al. [22] and Nagata et al. [36] described the dense seismograph array 
network installed at Chiba experiment station, Institute of Industrial Science, 
University of Tokyo in 1982. It measured earthquake records in three dimensions, 
which can be used in evaluating the spatial correlation of earthquake ground 
motions and in the study of soil amplification. Katayama et al. [24] had given a 
detailed description of Chiba array strong motion database and its effective use. 

Yamazaki et al. [49] provided information about location, site condition at Chiba 
experiment station, corrections for orientation errors, process of development & 
organization of database and emphasized the use of event codes in determining 
engineering properties of strong ground motion. 

Tsujihara et al. [47] studied the subsurface soil properties based on Chiba downhole 
array data and formulated a ‘frequency domain problem’ and implemented the ‘transfer 
function method’ in obtaining the shear wave velocity (Vs) & quality factor (Q). Estimated 
shear wave velocity (Vs) was less than the elastic wave velocity (Fig. 3.2). Events with 
varying PGA were analyzed and the study reported that the soils below 10m depth are 
relatively stiff (e.g. Vs > 300 m/sec) and no non-linearity was observed even during large 
excitation of the ground. The quality factor (Q) obtained was different in every event, so that 
the identification could hot be carried out accurately to make a definite conclusion. 
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Figure 3.2: Profiles of identified Shear wave velocity and Quality factor 
(Tsujihara et al. [47]) 

Sawada et al. [44] formulated multiple reflection theoiy in frequency domain to estimate the 
dynamic soil properties. Chiba array records were analyzed and identified Shear wave 
velocity (Vs) was in reasonable agreement with that of PS logging test (Fig. 3.3). However, 
there was a large scatter in identified Quality factor (Q) values indicating the less reliability 
of the employed method. 
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Figure 3.3 : Estimated values of Shear wave velocity and Quality factor 

(Sawada et al. [44]) 



Katayama et al. [28], in his Soil amplification study report on Chiba array 
downhole data, showed that the soft soil layer few meters (0 - 5 m) below the 
ground surface got amplified more than the rest of the stiff layers below it. 
Katayama and others [23, 25, 26, 27] used the three dimensional nature of Chiba 
array records to study the behavior of wave propagation, spatial correlation and 
affect of separation distance between two stations. 

Ghayamghamian et al. [17] used downhole ground motions recorded at five vertical sites in 
Japan to investigate the nonlinear site response and in-situ dynamic soil properties. The 
hysteresis behavior of soil was found to be changed for different time windows (i.e. for 
different levels of shaking) (Figure. 3.4) and determined nonlinear behavior of the site. The 
shear moduli inferred from the events producing strong ground motions were substantially 
lower than that producing weak ground motions. 
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Figure 3.4: Shear stress verses Shear strain for different time windows 
(Ghayamghamian [17]) 








Two records of Chiba vertical array where PGA recorded was 33.5 gal and 301.1 gal were 
used. Shear stress-strain hysteresis loops were plotted for earthquakes with PGA 301.1 gal (in 
9-14, 14-19 & 19-24 sec time windows) and PGA 33.5 gal (in 10-15, 15-20 & 20-25 sec time 
windows) and shown in Figure 3.4. A mild detectable non-linear behavior was recorded 
during the event with PGA 301.1 gal, which could be observed by changing shear modulus 
and damping ratios in different time windows. However, this phenomenon was not observed 
during PGA 33.5 gal due to the small levels of shaking. 

Katayama et al. [28, 29] investigated the soil amplification characteristics in terms of peak 
ground acceleration and transfer function based on the Chiba array observation records. The 
recorded peak accelerations at borehole CO at GL -1, -5, -10 and -20 m were normahzed by 
those at GL -40 m for the 27 events (Chiba database) to study the amplification of peak 
acceleration. The ratios for the three components i.e. East-West (EW), North-South (NS) and 
Up-Down (UD) for individual events are shown in Figure 3.5. 
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Figure 3.5: Amplification ratios of peak ground acceleration for 27 events at 
borehole CO (Katayama et al. [29]) 







Amplification took place mostly with in the few meters below ground level due to top 5m soil 
and the underlying layers are relatively denser and not subjected to considerable 
amplification. Transfer functions of two events (8722 with PGA 327.1 cm/sec^ and 8725 with 
PGA 24 cm/sec^) were compared to study the possible effects of non-linear response of soil 
deposits on the transfer function. Smoothened transfer function for these two events is shown 
in Figure 3.6. Effect of non-linear response of soil was observed to be small for the soil 
below 10 m depth due to the denser sub soil layer. 






Figure 3.6: Smoothened transfer functions estimated from NS component of the 1987 
east off Chiba prefecture earthquake and one of its after shocks in 
borehole CO (Katayama et al. [29]) 







3.5 Other studies on Downhole arrays 


Yeon - Tein et al. [51] described the strong motion array in Haulien, Taiwan called SMART- 
2 and development of LSST (Large Scale Seismic Testing) experiment on hard site in the 
Haulien area. SMART-2 was designed for studying the rupture process of earthquake faulting 
and near field ground motions. Haulien LSST experiment was designed for soil structure 
interaction research. 

Ansary et al. [3] discussed the characteristics of micro tremors based on the spectral 
analysis on Haulien LSST microtremor data. The amplification characteristics were 
observed in the free field as well as in the back fill. The amplitude ratio of theoretical 
rayleigh wave matching with the observed amplitude ratio of microtremors clearly 
indicated the influence of rayleigh wave except for some low and high frequency 
microtremors which may consists of low waves. Back fill soil microtremor records did 
not reflect good matching of amplitude ratio between theoretical rayliegh wave and 
observed, alluding the difficulty in predicting the characteristics of micro tremors in the 
case of back fill soil. 

Cheng et al. [12] used Haulien vertical firee-field arrays to study the response of soil during 
strong motion earthquake. The geological profile is very typical consists of loose sands for 
the top 5 ni underlying gravel deposits of unconsolidated conglomerate composed of pebbles 
varying in diameter 10 cm to 20 cm. The assumption was that the horizontal ground motions 
were because of vertically propagating Shear wave and hence One-dimensional shear beam 
model was employed. 

Amplitude spectrum and phase spectrums were estimated from spectral analysis to calculate 
the modal frequency from amplitude spectrum and phase angle from phase spectrum. Single 
and two layer characteristic equations for transfer function were derived from wave 
propagation theory been used to calculate the parameters for system identification by 
substituting the predominant frequencies obtained from amplitude spectrum. The velocities 
obtained in EW direction between different layers were lower than the velocities obtained in 
the NS direction, clearly indicated the anisotropic behavior of the soil. 



Archuleta et al. [4] conducted site amplification study on Gamer valley downhole array and 
observed the modification in the amplitude of the acceleration in the surface layers, which is 
basically soil overlain by weathered granite. Gamer valley downhole array of accelerometers 
was used to measure the effects of local site conditions on the amplification and attenuation 
of seismic waves. 

Archuleta et al. [5] studied the Gamer valley accelerometer records for amplification 
phenomenon. The amplitude of waves will not be effected by the layer when the thickness is 
less than the wavelength, moreover the amplification occurs when the frequencies whose 
quarter wavelength is less than the layer thickness. Analysis shown that larger earthquake 
was amplified less than the smaller earthquakes. This is due to the lower comer frequency in 
the larger earthquakes, as the comer frequency decreases, the wavelength increase and hence 
the amplitude decreases. 

Pecker et al. [38] used the weak motion records of Gamer valley array to validate the low 
strain soil properties obtained from soil survey at the experiment site. Spectral ratio technique 
was used to analyze the weak ground motion records, consequently to compare with the 
results obtained from elastic behavior from soil survey. The soil behavior was isotropic and 
yielded in reasonable validation of results from the geotechnical survey. 

Peter et al. [40] analyzed seismic signals using multiple signal characterization (MUSIC) 
method. Results showed that the method was capable of resolving multiple closely spaced 
source and work well with both stationary and non-stationary signals. 

Yasuyuki et al. [50] investigated- the behavior of embankment model resting on sandy 
ground, using the measured records from buried instruments in the model. Pore water 
pressure was found to be proportionate to the input acceleration characteristics. Ground 
beneath the embankment has not undergone cyclic mobility effect unlike it was found in free 
ground and large settlement was observed in the embankment due to vibrational earthquake 
motion than a shake-t 5 q)e motion. 



Toshimi et al. [46] defined engineering bedrock as a linear behavior regime even during the 
strongest shaking after studying the bedrock waves from the Sendai region borehole 
observation network, Japan, using 1-D soil modeling. Such waves estimated from site 
amplification were observed to be much lower than the observed surface & borehole records 
themselves. Therefore they could be used in aseismic analysis of important structures and 
used in the study of source characteristics of the observed events as well as the path effects 
between the seismological bedrock and the engineering bedrock. 

Safak [42] has given discrete time wave propagation analysis to investigate the site 
amplification in layered media. The model uses only three parameters for determining the site 
amplification, viz. R (reflection coefficient), r (one way travel time of waves in the layer), 
and Q (quality factor). This method is much simpler and more accurate than the frequency 
domain techniques and useful in applications such as random vibration studies, site effects on 
response spectra and in identification of site characteristics fi'om recorded motions. 


3.6 Summary 

From the review of literature carried out directly related to the current problem, it can be 

summarized that, 

• Tangible and reasonable information on dynamic site behavior can be obtained from the 
study of downhole array records. 

• Downhole array observation is a useful tool in estimating the in-situ dynamic properties. 

• Direct analysis of the downhole records can overcome the constraints involved in 
preparation and testing of laboratory procedures to simulate the exact field conditions. 

• Downhole array earthquake records are powerful tools that complement lab testing to 
study the dynamic behavior of soil. 



Chapter 4 

CORRELATION ANALYSES 


4.1 Introduction 

Correlation analyses have several useful applications in the area of seismic wave propagation 
such as, 

• Identifying similarities in signals, 

• Determination of travel times (time lag), 

• Identifying replicas in noisy signals and 

• To characterize multiple travel paths, lossy and dispersive media. 

Several other researchers had earlier employed Cross-correlation techmque, to estimate the 
time lag. Few among them that directly related the current study are; 

• Campanella et al.[8] in seismic cone penetration testing using frequency domain 
method 

• Elgamal et al. [14] in estimating site dynamic properties 

• Matsunaga et al. [35] in calculating phase velocities. 


As was indicated in the review of literature, in identifying sub surface shear wave profile 
(using time & frequency domain techniques), Cross-correlation technique is a very reliable 
procedure and enjoyed popularity among the researchers (Santamarina et al. [43]). 

In this study. Cross-correlation technique in time domain method is used in estimating the 
average Shear and Primary wave velocity distribution profiles between the downhole stations 
at different depths. The advantage of this technique is its ability to identify the replicas of a 
signal in other signals, even in the presence of significant background noise. These profiles 
were then compared with those obtained from the Geo-physical measurements (Figure.2.3). 



4.2 Methodology 


The purpose of employing Cross-correlation in the present work is to identify the time 
taken by the wave to travel between two downhole stations. If a.i and Uj are two 
acceleration histories recorded at two different stations i and j, then the Cross-correlation 
in continuous time is a function of the time shift r and is given as: 

00 

= \aXt).aj{t + T).dt ( 4 . 1 ) 

-CO 

where, a/t) and a/t) represent the time records of wave passing through the first and second 
receivers placed z m apart (in depth) as shown in Figure. 4.1. r represents the time delay 
between the two records and t is the domain of integration. However, in reality, limits of 
integration in the time domain will be finite and above Cross-correlation can be normalized, 
to keep the maximum absolute value of Q to be less than one. 



a, ft) 


a,- (t) 


Figure 4.1: The recorded signals at z m apart for time domain Cross-correlation 

Cross-correlation fimction between downhole acceleration time histories u, (0 and {t) 

recorded at stations i and j may be expressed as: 
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N-m^i 


I 1 A' 1 JV 

— Yuf(MAr) — ^af(nAt) 


(4.2) 


Here, time delay, r = m-At, where m = 0,1, 2,3 N. N represents the total number of 

records and At is the time step of digital data. 



In principle, the Cross-correlation of signals at adjacent depths is determined by shifting the 
lower signal, relative to the upper signal, in steps equal to the time interval between the 
digitized points of the signals (Campanella et al. [8]). 

At each shift, the sum of the products of the signal amplitudes at each interval gives the 
cross-correlation for that shift. After shifting through all the time intervals. Cross-correlation 
function verses Time-shift can be plotted, and the time shift interval is used to calculate the 
interval velocity. 

The process is shown schematically in Figure 4.2, where the lower signal has been shifted to 
the left giving the maximum correlation. 



Figure 4.2: Cross-correlation method to obtain Time interval ( ) 


Cross-correlation function can be calculated in frequency or time domain. In the present 
work, ‘time domain approach’ is adopted and the process of Cross-correlation calculation is 
shown in Figure 4.3. 
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Figure 4.3: Cross-correlation calculation in time domain 


Cross-correlation function reaches a major peak at a time delay T=Td, where Td is the time for 
seismic waves to travel from station i to station j. Then, the apparent velocity Vd of wave 
propagation between stations i and j may be estimated as, 


Vd 


d 


(4.3) 


where, d is the known distance between i and j stations. 









Acceleration frequency content records for each event can greatly influence the shape of the 
Cross-correlation function. Extraneous noise at the input or output reduces the relative 
contribution of individual correlation peaks and increases the random error in the analysis, 
but it does not distort or bias the results. 


Cross-correlation function obtains sharp peaks and is effective, when the low frequency 
content (f < 1/ (2 Xd)) were filtered out, to distinguish the incident & reflected peaks. Such 
filtering was necessary to allow separation of incident and reflected wave peaks, in view of: 


• the relatively narrow bandwidth of large amplitude frequency components and 

• the relatively close spacing of upper most accelerometers. 


Assuming the time required for the' wave to travel from seismograph u,- to surface is ^d(a,-*s) > 

where a, refers seismograph station at the specified depth, i refers to the station index from 
surface (1,2,3 . . .n stations). The velocity between the station to the surface is given by: 


V. 


s(a, s) 


-* s) 
^ d(a, -¥ s) 


(4.4) 


where, is the depth of the station from the surface and is the velocity 

between the station to the surface. Therefore the average velocity between any two 
consecutive stations is given by: 

«''>«•«. (4-5) 
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(4.6 a) 
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4.3 Characteristics of Cross-correlation Function 

• The value of shift for the peak (positive or negative) in Cross-correlation indicates 
the time delay between the replicas. 

• If the maximum Cross-correlation function has a higher value, the peaks in both 
the records have similar magnitudes and if it has a smaller value, the two peaks 
have large variation in magnitude and their wave shapes are different. 

• When the Cross-correlation function is maximum, time lag is represented by the 
phase difference between the two records. 

• The peak close to the origin represents the case where transducer is very close to 
the discontinuity such as free surface or a close spacing between transducers. 

• Cross-correlation function magnitude represents the energy carried by the waves. 
Reflected waves have smaller magnitude, as some energy is lost due to damping. 

• If the signal and its replica were of opposite polarity (-a,), the maximum value of 
Q ^ would be negative. Polarity reversal can be caused by the wiring of transducers, 

or physical reasons such as the reflection of a signal at a free boundary. 

4.4 Limitations 

Cross-correlation technique provides average shear wave velocity between two downhole 
stations. In general, the estimates of shear wave velocity are reliable only if, 

• The earthquake records are of broad frequency bandwidth, 

• Propagation velocities allow for sufficient separation between incident and reflected 
wave components, 

• Low amplitude excitation prevails so as to avoid non-linear soil response and 

• Signal is dominated by body shear waves (rather than the surface waves). 


44 



4.5 Results and Discussions 


Cross-correlation technique is reliable when two stations recording acceleration history are 
reasonably apart. If they are close by, it is very difficult to judge the time lag, because the two 
records may be similar with in the short distance. In addition, data recorded close to the 
surface might be effected by the domination of the surface waves, resulting in poor 
distribution of Cross-correlation function in estimating time lag in 0 - 5 m layer. Therefore, in 
the present study, time lag was estimated between 0-5 m, 0-10 m, 0-20 m and 0 - 40 m 
layers rather than in 0 - 5 m, 5 - 10 m, 10 - 20 m, 20 - 40 m layers. These results were then 
used to calculate the time taken by the waves to travel in 40 - 20 m, 20 - 10 m, 10 - 5 m and 
5 - 0 m layers. 

Normalized Cross-correlation function is calculated using Equation 4.2. Time lag (r^) is 
calculated fi-om the lagged cross-correlation peak from the center. The firequency components 
( / < l/(2r^) ) were filtered-out fi’om the signal 1 and signal 2 to clearly identify the incident 

and reflected wave peaks. Cut-off fi'equency values for the records used in this analysis are 
shown in the Table 4.2 & 4.3. From the Cross-correlation diagrams (Figure 4.4 - 4.13), it is 
evident that the time lag is decreasing with the depth and vice versa, i.e. the Shear wave 
velocity decreases as the wave reaches the surface. It can be observed that some records were 
not used in the calculation of Cross-correlation function especially between the surface and 5 
m layer. The reason being this layer is dominated by the surface waves rather than the body 
waves. 

Spectral density function between two time history records may be defined as the Fourier 
transformation of the correlation function between those records. Good frequency bandwidth 
distribution can result in getting accurate estimation of cross-correlation function and the time 
lag. Figures 4.14 - 4.16 illustrate the Cross-spectral density function, in which, the 
earthquake records did not contain enough firequency bandwidth for a meaningful 
interpretation of Cross-correlation function, (e.g. the earthquake 8823 for which very narrow 
fi'equency bandwidth (3-4 Hz) could be observed almost at all depths especially in the NS 
direction). Such narrow bandwidth of recorded signal produced an unreliable Cross- 
correlation function and was not used in the analysis. 
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Figures 4.17 - 4.19 shows the Cross-correlation function used to calculate Primary wave 
velocity. The amplitude of the function increased with depth and two main peaks can be 
observed as the incident and reflected waves. As can be seen from the Figures 4.4 — 4.13 and 
4. 1 7 - 4. 1 9, that the time lag increased with the depth. The time lag values obtained from the 
cross-correlation analysis for different depths to the surface in the calculation of shear wave 
and primary wave velocities are shown in Table 4.4 and 4.5. The calculated shear and 
primary wave velocities between the consecutive stations are shown in Table 4.6 and 4.7 
respectively. Shear and Primary wave velocity profiles were estimated for different events by 
Cross-correlation analysis and are compared with geophysical measurement velocity profile 
(Figures 4.20 - 4.21). These velocities in EW & NS directions are in reasonable agreement 
with the velocities measured using in-situ geophysical method. 

4.6 Summary 

• Cross-correlation method is a reliable technique in estimating time lag when the 
recording stations are far apart. 

• Filtering was effective in clearly differentiating the incident and reflected wave peaks. 

• Estimation of wave propagation velocities were reliable when the signal is dominated 
by body waves and having large frequency band width. 

• The estimated shear wave velocities for the events where the PGA < 41 cm/sec^ (low 
amplitude excitation) recorded on the surface, were reasonably in good agreement 
with the geophysical measurements. 

• Shear wave velocities were observed to be higher than the geophysical measurements, 
especially in 10-20 m and 20-40 m layers. 

• Some of the events were not employed to estimate shear wave velocities due to 
narrow frequency bandwidth. 
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The estimated shear wave velocity during an event where recorded PGA (Peak event 
8722, PGA is 327.1 cm/sec^ on the surface) is observed to be in agreement with the 
geophysical measurements.' However, for such a large PGA, significant non-linear 
response is expected and shear wave velocity should be lower. As observed by other 
researchers and corroborated in this study, non-linear response of soil is small even 
during peak PGA due to rigid sub soil layers and hence the velocities were almost 
coinciding in both higher and lower PGA studies. 



Table 4.1: Lower cut-off frequencies employed in the calculation of cross-correlation 

(Shear wave velocity estimation) 
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Table 4.2; Lower cut-off frequencies employed in the calculation of cross-correlation 

(Primary wave velocity estimation) 



8903 13.2 









































































































Cross-correlation 





Time Lag (sec) Time Lag (sec) 


Figure 4.5: Cross-correlation of accelerafion iisfories during the earthquake 8510 
at borehole CO 
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Figure 4.10: Gross-correlation of acceleration histories during the earthquake 8722 
at borehole GO 
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Figure 4.14: Cross-spectral density function during tie eartlquake 8307 
at borehole CO 
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Figure 4.15: Cross^spectratdensity finictioir during the earth^iiake 8406' 

at borehole eo 















Cross spectral density amplitude 


EW Direction 





01 23456789 10 


01 23456789 10 


01 23456789 10 



Frequncy (Hz) 


Figure 4.16: 
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Figure 4.17; Cross-correlation of acceleration histories during the earthquakes 
8722 and 8723 at borehole CO 















































Table 4.5: Estimated shear wave velocities between the borehole stations 
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Table 4.6: Estimated Primary wave velocities between the borehole stations 
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Figure 4.20: Estimated Shear wave velocity using Cross-correlation method 
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Figure 4.21: Estimated Primary wave velocity using Cross-correlation method 
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Chapter 5 

SHEAR STRESS - STRAIN HISTORIES 


5.1 Introduction 

Calculating site shear stress-strain histories directly from the recorded acceleration histories 
enables to; 


• understand the site amplification (linear or non-linear) under earthquake loading 

• estimate the stiffness of the soil 

• estimate the hysteretic damping of the soil 

• estimate the induced strain in the soil. 

Since the loading amplitude on soil varies with time during the dynamic loading, soil 
stiffness changes for each cycle. It influences the wave propagation velocity through soil 
media and is inversely proportional to the induced strain in the soil. Ground strain is effected 
by ground acceleration amplitude and ground displacements. In the present study One- 
dimensional shear beam idealization is used to describe the site seismic lateral response in 
line with other researchers. Similar studies by other researchers include: 

• Abdel - Ghaffer et al. [2] employed a single degree of freedom system to calculate 
the shear moduli and damping values of Santa Felicia earth dam. 

• Koya & Matsao [30] employed a force-relative displacement relation to study the 
results of shake table tests. 

• Zeghal et al. [54] used One-dimensional shear beam idealization to study Lotung 
downhole array soil non-linear properties by calculating shear stress-strain 
histories. 

In the current study as well as those reported in the literature review, longer period waves 
(< 0.25 - 0.35 HZ) were filtered out in calculating the ground stress-strain to estimate the soil 
d 5 aiamic properties. 



5.2 Methodology 


Using a shear beam model to describe site seismic lateral response, shear stress t (z, t) at any 
level 2 can be expressed as: 

rfeO = I pu (CAdi (5.1) 

0 

where, /?= mass density, and w (j^J) - absolute acceleration at level C 

Using the linear interpolation between downhole accelerations, the discrete counterpart of 
shear stress at level 2 is shown in Figure 5.1. 



Figure 5.1: Site Discretization 

From the Newton’s Second Law of Motion (F = ma) and assuming unit area of shear beam in 
the calculation of stress, 
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From the above equation the stress at any level z, is expressed as, 
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Like -wise the force at the mid point of levels 1 and 2 is given as. 
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Substituting u j from equation 5.5, 
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Stress is given as, r , (r) 
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From the above equation the strain at any level (z,.i + Zi)/2 is expressed as. 
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In equations 5.4 and 5.9, i = 2,3 . . . and i and i-1/2 refer to levels z , (of the i* accelerometer) 
and (z,.i + Zi )/2 (halfway between accelerometers i and (i-1)) respectively, t, (t) = t ( z,, t), 
z'q = Az, is spacing interval between accelerometers. 



These stress estimates are second order accurate. And the corresponding second-order 
accurate shear strain between the levels 1 and 2 levels may be expressed as, 
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From the above equation the strain at any level (zi-i + z,)/2 is expressed as, 
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Similarly strain at level 2 may be expressed as. 
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From the above equation the strain at any level z, is expressed as, 
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In equations 5.11 and 5.13, i = 2,3... and u, = u (z„ t) is absolute displacement evaluated 
through double integration of the recorded acceleration history u (z,. ,t). 



5.3 Dynamic Site Properties 


The evaluated shear stress strain histories at Chiba were used to estimate soil stiffness and 
material damping using the model, 

r - G/+7J/ (5.14) 

in which G and p are shear modulus and viscous damping. 


These stiffness and damping parameters during a given time window (f, tf) may be evaluated 
using a least squares fit to minimize the difference between shear stresses estimated from 
acceleration histories (Equations 5.15 and 5.16), and those predicted by linear model of 
Equation 5.14. 


Thus, least squares minimization dictates: 






If If 


^Tydt jyydt - jrydt jy/dt 


0 f ‘j Y 

jyydt jyydt - jyydt 


( 5 . 15 ) 


(^1 5 ^/ ) 


I, t, 




^rydt ^yydt ~ ^tydt ^yydt 


ft, 


^yydt ^yydt — ^yydt 


J 


(5.16) 


where x is stress history, y is strain history and /is strain rate history. 
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5.4 Results and Discussions 


Site shear stress-strain responses were evaluated at four locations within the instrumented 
zone. Shear stress-strain hysterisis loops were shown in Figure 5.14 - 5.19. Overall, the 
employed linear model is seen to capture the essential features of recorded response. 
Calculated shear wave velocities and damping ratio values at different depths for the 25 
events are shown in Table 5.1. Shear modulus & damping values were calculated using the 
Equations 5.15 and 5.16 and shear wave velocity was calculated based on the shear modulus 
mass density. 

Displacement histories were calculated by double integrating the filtered acceleration 
histories and used to estimate the strain histories. These strain histories were calculated using 
the Equations 5.4 and 5.9 at halfway between accelerometers at 2.5, 7.5, 15.0 & 30.0 m 
depths in EW & NS directions and shown in Figure 5.2 & 5.9. These figures clearly show 
that shear stress was increasing with depth and the shear strain was decreasing with depth. 
Ground strain at surface (0 - 5 m layer) is significantly large increased compared to other 
layers, clearly indicating the surface soil amplification during the earthquake loading. 

The optimized stresses calculated by assumed linear model with damping (Equation 5.14) are 
compared with the shear stress calculated using the one-dimensional shear beam model 
(Equations 5.4 and 5.9). The optimized shear stresses are in good agreement with the 
recorded shear stresses. These figures also show the variation in the stiffness of the material 
at various depths below the ground surface. The surface layer of loam is subjected to more 
ground strain compared to silty-clay and sand at other depths. 

The calculated shear wave velocities are compared with the geophysical tests as shown in 
Figure 5.20. Shear wave velocities are found to be slightly lower than those obtained by 
geophysical method due to reduction in soil stiffness and mild non-linear behavior. Damping 
value for the top layer was also higher compared to other depths (Table 5.1). This 
phenomenon is supported by the Figures. 5.2 - 5.9. From these figures, It can be observed 
that shear stress was lowest at the top layer and increased with depth where as the shear strain 
was maximum at the surface and decreased with depth. 
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Summary 


The slight difference in 'the estimated shear wave velocities & geophysical 
measurements indicate that geophysical measurements induced much lower strains 
compared to earthquakes used in the analysis. 

Shear stress levels had increased with depth but the shear strain decreased with depth. 

Amplification in the top layer had produced high excitation levels and generated high 
strains due to loss of stiffness, which is supported by lower bound values estimated by 
the stress-strain imaging. 

The ranges of values of shear wave velocity profile estimated using this method were 
in good agreement with geophysical measurements except for top 5 m layer. It further 
corroborates that the behavior of soil at Chiba during the earthquake analyses was 
mildly nonlinear even during the strongest ground motion. 



PGA Shear wave velocity Vj, (m/sec) and Damping ratio Tj/G (%) 

_ ^ ( cm/sec^ ) — — — 
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Figure 5.2; Computed stress histories at borehole CO during the earthquake 8510 
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Figure 5.3; Computed strain histories at borehole CO during the earthquake 8510 













Figure 5.4: Computed stress histories at borehole CO during the earthquake 8525 
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Figure 5.8; Computed stress histories at borehole CO during the earthquake 8726 












Figure 5-9: Computed strain histories at borehole CO during the earthquake 8726 
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Figure 5.10; Evaluated and optimal stress histories at borehole CO during the 
earthquake 8510 
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Figure 5.1 1 : Evaluated and optimal stress histories at borehole CO during the 
earthquake 8525 
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Figure 5.12; Evaluated and optimal stress histories at borehole CO during the 
earthquake 8722 
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Figure 5.13: Evaluated and optimal stress histories at borehole CO during the 
earthquake 8726 
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Figure 5.14: Computed shear stress-strain histories at borehole CO during the 
earthquake 8510 
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Figure 5.17: Computed shear stress-strain histories at borehole CO during the 
earthquake 8726 










Shear stress (k Pa) 




Figure 5.18; Computed shear stress-i 
earthquake 8806 
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Chapter 6 

TWO-LAYER SOLUTION 


6.1 Introduction 

Strong ground motions at certain frequencies are largely affected by the amplification 
characteristics (shear wave profile) of subsurface layers. Amplification characteristics can 
provide information for the prediction of ground motion for a specific site and for the seismic 
safety assessment of existing structures. 

In the present study, spectral ratio technique was used to study the response of the sub surface 
soil layers at borehole CO of Chiba array. Two-layer solution was used to estimate the 
fundamental frequencies of soil layers and shear wave velocities between consecutive 
borehole stations. 

Some other studies where spectral ratio techniques were employed to study amplification 
characteristics are: 


• Chang et al. [9] in ground response analysis at Lotung downhole array. 

• Chang et al. [10] in back calculating equivalent shear moduli (and their variations 
with shearing strains) from the resonant frequency of the Fourier spectral ratio at 
the Lotung site. 

• Katayama et al. [28] in studying the ground amplification and calculation of 
natural frequencies. 

• Wen et al. [48] studied the soil nonlinear seismic response by comparing the 
spectral ratios of surface to downhole horizontal accelerations on weak and strong 
motion based on the array records of SMART land SMART 2 in Taiwan. 

• Huey-chu et al. [21] to investigate the site amplification of four level Dahan 
downhole recordings. 
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6.2 Methodology 


The transfer function between the acceleration histories recorded at two depths can be 
calculated to investigate the characteristics of ground vibrations and to identify the modal 
frequencies of vibration of the soil layers. The transfer function was calculated between 
surface and other depths. 

If a, and aj are surface and the depth at which ground motion is recorded (Figure 6.1), then 
Fourier spectral ratio is given as: 


H 


a,aj 


if) 


H., if) 
H.^if) 


( 6 . 1 ) 


where (/) and are the Fourier transforms of the ground motion a, and aj 



Figure 6.1: Recorded signals at z m apart for the calculation of transfer function 


//„„ (/) represents the transfer function of frequency / and ^ (/)]j defines the 

amplitude spectrum. Arc tangent value of the ratio between the imaginary and real parts of 
the transfer function defines the phase spectrum. Amplitude and phase spectrums are used to 
identify the natural frequencies of vibrations of soil layers. 

Since the spectral ratios calculated from the earthquake records are random in nature, 
Hanning window technique was used to smoothen the amplitude function and to locate the 
peak values for modal frequency identification. 
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The averaging operation is done on the spectral ratio function to calculate the smoothened 
function y (/) is given as: 

y,(/) = 0.25y,.,(/)+0.5y,(/)+0.25y„,(/) (6.2) 


where, i vary from 2 to N-1. 

Smoothening procedure is repeated for the required number of times (8 to 10 times). The 
lowest frequencies at which the Fourier spectral ratio reaches its peak and where phase angle 
is 90° are identified. Equivalent / average shear wave velocities for each depth-interval were 
back calculated using resonant frequencies of the Fourier spectral ratios in conjunction with 
Two-layer solution derived by Madera (1970). 


The equation of motion for a vertically propagating shear wave is: 


P 


d^u 

W 


dz \ dz) 


where. 


(6.3) 


p = Mass density of soil 
X = x{z,t) = horizontal soil displacement 
z - vertical coordenate 
t = time 

G= pv^ ~ shear morulas of soil 
Vj = shear wave velocity 


The solution of the above Equation with boundary condition of zero relative displacement at 
the bottom boundary of the layers A and B, and zero shear stress at the ground surface is: 
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tan 


(6.4) 


U/J 


r 

tan • 

V 


r \ 

!lL 

2/J 


Pb ^ b fs 
Pa H, L 


Where ^Fundamental frequency of layers A and B respectively 

f = Fundamental frequency of the total soil layer 
PaiPb ~ densities of layers A and B respectively 
H^,Hg = Thickness of layers A and B respectively 


With the known fundamental resonant frequencies of layer A and the two layer profile (/a and 
f respectively), fundamental frequency (ff) of layer B is computed from Equation 6.4. Shear 
wave velocities of layers A and B are computed by Vsa = 4Hyi^ and Vsb = 4Hq/b for a 
system of three layers of soil on top of an elastic half-space as shown is Figure 6.3. 

The two upper layers are regarded as an equivalent uniform single layer of thickness h^ and 
mass density where, 

K=K + hg (6.5) 

pA^i^hPA+KPs)! iK+h) ( 6 - 6 ) 


Upper 3 layers of soil are regarded as an equivalent two-layer system of thickness and he 
(or and ho...), mass density and p^ and resonant frequency (is fundamental 

frequency of the total soil layer above layer c) and^, where^ is calculated using Equation 
6.4. This procedure is followed for n number of layers and the velocity in each layer is 
calculated using 4HQ^ where Hn,^ are the thickness and fundamental frequencies of n* layer. 
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Figure 6.3: Models used for calculation of Shear wave velocity 
(a) Single layer system (b) l’wo-la\ cr system 
(c) Three-layer system (d) F(juivalent fw<»-layer system 



6.3 Results and Discussions 


Simplified ground model adopted for estimating the shear wave profile the current analysis is 
shown in Figure 6.4. 
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GL - 20.0 m 


-■ a4 


Gd Pd ho 

GL - 40.0 m ■ as 

ai, a 2 , as, 34 and as are the accelerometers at surface, 5 m, 10 m, 20 m and 40 m respectively 


Figure 6.4: Ground model used for system identification 


Transfer functions were calculated between surface and 5 m, surface and 10 m, surface and 
20 m and surface and 40 m for four events 8722 (PGA 327.1 cm/sec^), 8706 (PGA 14.0 
cm/sec^), 8903 (PGA 28.9 cm/sec^) and 8519 (PGA 82.2 cm/sec^) and were compared in 
Figure 6.5. 



As, csn be seen from Figure. 6.5, peaks in transfer function between suiface and 5 m and 
surface and 10 m exhibited some shifting with amplitude of PGA while this effect was not 
clear for other depths. This shows that there is no significant nonlinearity in the response for 
different input motions. 

Table 6.1 presents the identified modal frequencies. EW and NS components of amplitude 
and phase spectrums of selected events are shown in Figures 6.6 — 6.15. Significant soil non 
linearity was not detected in layers below 10 m even during the peak event (PGA 8722, 327.1 
cm/sec^) due to the denser soil layer and the strain in these layers was not so large to soften 
the soils. Identified shear wave velocities are consistent with the geophysical measurements 
as shown in Figure 6.26. 


6.4 Summary 


• Fundamental frequency depends on PGA and is affected by the response of soil deposit. 

• Fundamental frequency was higher in the top layers and decreased with depth. 

• Changes in shear modulus or shear wave velocity were implied by the variation in 
fundamental frequencies and amplitude of ground acceleration. 

• Reduction in fundamental frequency decreased with increasing thickness of the soil 
column. 

• Shear wave velocities estimated during the earthquake excitation for the top layer (0-5 
m) were substantially lower (2 - 37% variation) than the geophysical measurements 
implying the strains induced in the soil were higher than geophysical investigations. 

• During the events with strong ground motions (>0.35g), shear wave velocity profile was 
consistent with the geophysical measurements implying no detectable nonlinear behavior. 
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Table 6.1: Calculated Fundamental Frequencies and Estimated Shear Wave Velocities using Two Layer Solution 
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Figure 6.7; Amplitude spectrum and phase spectrum during the earthquake 8510 in 
north - south direction at borehole CO 
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Figure 6.12: Amplitude spectrum and phase spectrum during the earthquake 8717 in 
east - west direction at borehole CO 
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Figure 6.13; Amplitude spectrum and phase spectrum during the earthquake 8717 in 
north - south direction at borehole CO 
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Figure 6.16: Estimated Shear wave velocities by Two-Layer solution at borehole CO 
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Chapter 7 

CONCLUSIONS AND FUTURE RESEARCH 

7.1 Conclusions 

In the research reported herein, the recorded free field surface and down-hole acceleration 
histories at Chiba array were utilized to: 

• Evaluate shear wave propagation characteristics 

• Estimation of Shear wave and Primary velocity with depth 

• Evaluate site resonant frequencies 

• Evaluate site shear stress-strain histories 

• Estimate soil shear moduli and damping ratios 

Shear wave velocities were evaluated using correlation analyses and a two-layer solution 
technique. Resonant frequencies were estimated using spectral analyses. Identified shear wave 
profile was in good agreement with geophysical measurements. Evidence of soil softening 
during moderate to strongest earthquakes in the top 10 m of soil was also shown. 

A new non-parametric teclinique was used to directly evaluate soil seismic shear stress-strain 
histories from down-hole accelerations. The technique is based on interpolation between down- 
hole accelerations, and does not require the availability of bedrock or input motion records. The 
estimated shear stress and strain histories were employed to evaluate the soil shear stiffness and 
damping properties. The identified shear modulus and damping were able to capture the essential 
features of recorded response. There was a large scatter in identified damping in the top 10 m. 
This large scatter might be an actual reflection of the soft nature of top soil, along with other 
surface wave propagation characteristics that are not modeled by the employed simple one 
dimensional shear wave propagation concept. 
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A summary of the identified sliear wave velocities and compaiison with geophysical 
measurements is shown in Fig. 7.1. It is of interest to note that eaithquakc dciivcd shear wave 
velocities were biased toward the lower geophysical velocity data in the top 10 m. This might 
reflect the fact that geophysical measurements were performed at lower strain amplitudes as 
compared to the employed seismic records. It is also observed that below 10 m depth soil 
softening was not observed even during the event 8722 where PGA exceeded 300 gals. Other 
researchers using different teclmiques also observed this fact. Figure 7.2 compares the identified 
shear wave profile with the results reported by Sawada et al. [44], which corroborates the above 
finding. This is attributable to the stiff nature of the soil below 5 m depth. 

A thorough review of the literature and the results of this study suggest that the amount of 
information that can be inferred from the down-hole seismic records is a function of the site 
geology. The choice of the site is a function of the kind of infomiation required i.e., whether the 
required information is on propagation paths; radiation pattern; magnitude and fault rupture; or 
nonlinear / liquefaction site behavior. The identified properties show the potential of down-hole 
records as a valuable source of information on dynamic site behavior. Simple tools such as 
correlation techniques can be utilized to identify the low-strain properties. In addition, down-hole 
acceleration records can be utilized to evaluate the seismic shear stress strain histories, to analyze 
tire site response, and to estimate the soil dynamic properties such as shear modulus and damping 
properties. Finally it is emphasized that the techniques that were used herein may also be 
conveniently employed in conjunction with well-instrumented soil systems such as earth dams, 
retaining walls etc. 

7.2 Future Work 

The information and site properties provided by this study will be utilized in a irarametric system 
identification of Chiba site seismic response characteristics. The estimation will account for the 
two dimensional horizontal shear response. Parameter identification will be formulated on the 
basis of optimality criteria that related the observed and predicted seismic histories. Further 
research is required to fully quantify the effect of dense lower layers on site amplification. 
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Figure 7.1; Range of shear wave velocities estimated by various methods 
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Figure 7.2: Comparison of shear wave velocities with the results reported by 

Sawada et al, [44] 
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